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controlled  by  the  inviscid  equations.  The  equations,  written 
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are  approximated  by  MacCormack 1 s  second-order  accurate 
predictor-corrector  algorithm.  The  flow  tangency  conditions 
at  the  body  surface  are  satisfied  by  Abbett's  scheme  and  the 
outer  bow-shock  position  by  the  Rankine-Hugoniot  jump 
relations.  Internal  shock  waves  or  tangentially  discon¬ 
tinuities  are  captured. 

Comparisons  have  been  made  between  the  present  theory  and 
experimental  data  for  2.0  and  3.0  caliber  tangent  ogive  nose 
lengths  at  a  =  15°  and  20°  and  =  3.0.  Good  agreement  is 
obtained.  Results  were  also  obtained  for  a  5.0  caliber 
tangent  ogive  nose  length  configuration  for  comparison  with 
a  numerical  solution  obtained  by  solving  the  Navier-Stokes 
equation  and  with  experimental  data. 


Unclassified _ 

SECURITY  CLASSIFICATION  OF  THIS  PACEr*Fh«n  Dmtm  Entirmd) 

I 

it 


SUMMARY 


This  paper  addresses  the  problem  of  computing  the  steady 
inviscid  supersonic  flows  about  tangent  ogive  cylinders  at  high 
angles  of  attack.  For  such  cases  the  flow  on  the  leeward  side 
tends  to  separate  and  spiral  into  symmetric  vortices.  At  higher 
angles  of  attack  asymmetric  vortex  separation  can  occur.  In  this 
report  only  the  symmetric  cases  are  considered. 

Previous  analytical  studies  to  solve  these  flow  fields  have 
used  the  vortex-cloud  method  of  Mendenhall.  This  method  is 
basically  inviscid  and  semi-empirical.  Viscous  effects  have 
been  included  by  Pulliam  in  a  numerical  study  of  transonic  flows 
over  hemisphere  cylinders  and  supersonic  flows  over  tangent 
ogive  cylinders.  These  studies  were  based  on  the  Navier-Stokes 
equations,  the  former  on  the  time-dependent  three  dimensional 
Navier-Stokes  equations  and  the  latter  on  the  parabolized  Navier 
Stokes  equations.  To  obtain  a  more  efficient  procedure  we  use 
the  steady  Euler  equations  as  the  basic  governing  equations  as 
opposed  to  the  more  complicated  Navier-Stokes  equation.  The 
viscous  effects,  important  near  the  separation  lines,  are 
simulated  by  a  Kutta  condition.  The  separation  lines  are  deter¬ 
mined  from  experimental  data.  The  rest  of  the  flow  field  is 
essentially  controlled  by  the  inviscid  equations.  The  equations 
are  written  in  conservation  form  in  generalized  curvilinear 
coordinates.  The  equations  are  approximated  by  MacCormack's 
second-order  accurate  predictor-corrector  algorithm.  The  flow 
tangency  conditions  at  the  body  surface  are  satisfied  by  Abbett' 
scheme  and  the  outer  bow-shock  position  by  the  Rankine-Hugoniot 
jump  relations.  Any  internal  shock  waves  or  tangentially  dis¬ 
continuities  are  captures  by  the  scheme. 

Comparisons  are  made  between  the  results  based  on  the 
present  procedure  and  experimental  data  for  2.0  and  3.0  caliber 
tangent-ogive  noses  on  a  cylindrical  body  at  a  =  15°  and  20°  and 

1 


a  free-stream  Mach  number  =  3.0.  Good  agreement  is  obtained 
in  predicting  the  flow  field  and  the  vorticity  in  the  recircula¬ 
ting  region.  Results  are  also  obtained  for  a  5.0  caliber-long 
tangent  ogive  nose  configuration  for  comparison  with  a  numerical 
solution  based  on  the  parabolized  Navier-Stokes  equation  and 
with  experimental  data.  These  results  indicate  that  the  present 
procedure  gives  much  better  agreement  with  experimental  data 
than  the  Navier-Stokes  equation  results. 
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1.  INTRODUCTION 


As  a  result  of  past  work  (refs.  1  and  2)  done  at  Nielsen 
Engineering  &  Research,  Inc.  (NEAR)  for  the  Office  of  Naval 
Research  (ONR)  to  produce  engineering  prediction  methods  for 
missiles  up  to  angle  of  attack  of  50°  and  Mach  numbers  up  to 
3.0,  several  important  problems  in  high  angle  of  attack  aero¬ 
dynamics  have  emerged.  It  is  the  purpose  of  this  report  to 
describe  these  problems  and  to  attempt  their  solution  using  the 
Euler  equations.  While  these  problems  have  arisen  in  connection 
with  missile  aerodynamics,  they  are  equally  important  for  air¬ 
planes. 

The  present  work  is  a  logical  extension  of  the  work  NEAR 
has  carried  out  for  ONR  over  the  past  five  years.  During  the 
first  two  years,  paneling  (inviscid)  methods  together  with 
vortex  theory  were  used  to  develop  pressure  predictive  techniques 
for  complete  missiles  at  supersonic  speeds  and  angles  of  attack 
to  20°.  In  the  third  and  fourth  years  engineering  methods  were 
developed  to  predict  forces  and  moments  acting  on  canard  cruci¬ 
form  missiles  to  50°  angle  of  attack.  These  methods,  based  on 
data  base  and  rational  modeling,  produced  useful  methods  of 
reasonable  accuracy,  but  also  uncovered  important  aerodynamic 
effects  which  need  more  basic  study  in  their  own  right.  In  order 
to  extend  the  range  of  the  engineering  design  codes  and  to 
improve  the  rational  modeling  approach,  it  is  necessary  to  under¬ 
take  further  work.  It  is  believed  that  such  work  should  utilize 
recent  advances  in  viscous  and  inviscid  numerical  computation 
techniques  supplemented  by  further  experimental  work. 

Two  of  the  major  problems  which  have  been  uncovered  are  the 
adverse  effect  of  wing-body  interference  on  wing  lift  at  high 
angles  of  attack  and  the  special  behavior  of  body  vortices  at 
supercritical  crossflow  Mach  numbers.  The  former  problem  was 


7 


PRECEDING  page  BUNK -NOT  FIl*® 


covered  in  last  year's  contract  (ref.  3).  This  report  addresses 
the  latter  problem  only. 


We  will  first  discuss  the  background  of  vortex  shedding 
from  bodies  from  the  physical  point  of  view  so  that  the  signi¬ 
ficance  of  the  present  work  can  be  brought  into  focus.  Next 
the  state  of  the  analytical  art  is  reviewed  as  it  bears  on  the 
problem.  Then  the  means  of  introducing  compressibility  into  the 
problem  is  discussed.  This  means  is  the  same  as  that  used  for 
last  year's  contract,  namely,  solving  the  steady  Euler  equations 
of  motion  numerically.  The  viscous  effects  are  simulated  by 
imposing  a  Kutta  condition  where  the  flow  is  to  separate  from 
the  body.  The  numerical  procedure  will  not  be  described  here  as 
it  has  been  already  described  in  the  previous  year's  final 
report  (ref.  3) .  Changes  from  the  previous  method  will  be 
pointed  out  however.  The  report  concludes  with  a  section  on  a 
comparison  of  the  present  numerical  results  with  experimental 
data  (refs.  4,  5,  6,  and  7)  and  a  Navier-Stokes  solution  obtained 
by  an  implicit  numerical  scheme  developed  at  NASA/Ames  (ref.  18). 

2 .  BACKGROUND 

2.1  Vortex  Separation  for  Subcritical  Crossflow 

The  shedding  of  vortices  from  a  body  of  revolution  at  angle 
of  attack  has  a  number  of  regimes.  If  the  crossflow  Mach  number 
is  less  than  about  0.5,  the  crossflow  is  subcritical,  and  several 
types  of  steady  vortex  patterns  can  exist  on  the  leeward  side  of 
the  body  depending  on  its  fineness  ratio  and  its  angle  of  attack. 
The  two  types  of  vortex  shedding  are  characterized  as  symmetric 
and  asymmetric  as  shown  in  Figure  1.  The  boundary  between 
symmetric  and  asymmetric  vortex  shedding  for  subcritical  cross- 
flow  has  been  extensively  studied  by  a  number  of  investigators, 
and  a  good  summary  of  the  state  of  art  is  given  by  Chapman, 
Keener,  and  Malcolm  in  reference  9.  For  pointed  bodies  of 


revolution,  tnese  authors  show  that  asymmetry  starts  at  an  angle 
of  attack  of  about  twice  the  nose  semi-apex  angle.  For  a  conical 
nose  of  total  angle  25°,  asymmetry  might  occur  at  about  a  =  25°. 
For  an  L/D  =  3.5  tangent  ogive  nose,  asymmetry  starts  at  about 
32°  . 


2.2  Vortex  Separation  for  Supercritical  Crossflow 

If  the  crossflow  Mach  number  (M^sin  a)  is  much  greater  than 
about  0.5,  the  crossflow  is  said  to  be  supercritical.  In  this 
case  the  nature  of  the  vortex  separation  from  the  body  is  of  a 
totally  different  type  from  those  shown  in  Figure  1.  Consider 
Figure  2  in  which  the  shock  waves  are  not  all  shown.  Asymmetric 
vortices  do  not  now  appear  (ref.  9) .  Instead  two  elongated 
elliptic-like  regions  of  distributed  vorticity  form  on  each  side 
of  the  body  as  shown.  The  boundary  layer  separates  on  both  sides 
of  the  body  and  feeds  into  the  elongated  vorticity  regions.  As 
the  angle  of  attack  increases,  the  height  of  the  vorticity 
region  increases  as  shown  in  Figure  3  (ref.  2) . 

In  addition  to  these  features,  vapor  screen  photographs 
sometimes  show  what  appears  to  be  a  horizontal  shock  between  the 
regions  of  vorticity.  Also  there  appear  to  issue  from  between 
the  tops  of  the  regions  weak  alternating  vortices. 

2.3  Regions  of  Various  Vortex  Separation  Types 

Let  us  now  examine  the  a  -  diagram  in  Figure  4  which  shows 
the  regions  where  the  foregoing  three  classes  of  vortex  forma¬ 
tions  prevail.  If  Mc  =  0.5  is  the  boundary  between  the  sub- 
critical  and  supercritical  case,  then  above  the  boundary  we  have 
symmetric  elliptical  vorticity  regions.  Below  the  boundary,  we 
have  two  regions  separated  by  the  a  =  constant  line  between 
symmetric  and  asymmetric  separation.  What  this  figure  shows  is 
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that  the  supercritical  case  occupies  a  large  part  of  the  high 
angle  of  attack  region  over  the  Mach  number  range,  and  is  thus  of 
considerable  importance.  The  boundaries  as  given  in  this  figure 
are  based  on  approximate  rules  of  thumb,  but  the  general  con¬ 
clusion  is  not  changed  thereby.  Another  point  that  needs  to  be 
made  is  that  under  the  Mc  =  0.5  boundary,  the  symmetric  vortex 
model  is  significantly  influenced  by  compressibility  for  super¬ 
sonic  Mach  numbers.  It  is  the  effect  of  compressibility  on  the 
vortex  pattern  that  we  desire  to  predict,  particularly  that  for 
supercritical  crossflow.  No  method  for  studying  this  problem 
exists  except  possibly  computer  programs  based  on  the  unsteady 
three-dimensional  Navier-Stokes  or  the  parabolized  Navier-Stokes 
equations  (refs.  8  and  10) .  These  codes  require  very  long  com¬ 
puter  runs. 

2.4  Technical  Approach 

The  originally  envisaged  technical  approach  was  to  combine 
the  vortex-cloud  method  of  Mendenhall  (refs.  11  and  12)  with  the 
shock-capturing  Euler  code  described  in  reference  3.  The  vortex- 
cloud  methodology  has  presently  been  developed  for  incompressible 
flow,  and  the  flow  is  governed  by  the  Laplace  equation.  The 
vortices  are  introduced  into  the  flow  as  point  vortices  which 
also  satisfy  the  Laplace  equation.  Linear  superposition  is  used. 
In  the  case  of  compressibility,  the  general  flow  is  governed  by 
the  Euler  equations  which  permit  compressibility  and  shock  waves. 
The  Euler  equations  also  permit  distributed  vorticity  throughout 
the  field  as  long  as  the  flow  is  inviscid.  This  means  that  the 
vorticity  is  convected  but  not  diffused  by  molecular  action. 

This  assumption  should  be  quite  good.  We  now  discuss  the  vortex- 
cloud  method  for  incompressible  flow  and  then  describe  the 
necessary  modifications  and  complications  for  compressible  flow. 

2.4.1  Incompressible  Vortex-Cloud  Method.-  In  the  vortex- 
cloud  method  an  analogy  is  made  between  three-dimensional  steady 
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flow  about  a  body  and  two-dimensional  unsteady  flow  about  a 
circle.  Perhaps  a  good  way  to  envision  this  analogy  is  to  con¬ 
sider  plane  AA  in  Figure  2.  Let  the  plane  remain  fixed  and  let 
the  body  pierce  the  plane  at  time  t  =  0.  At  this  point  there  is 
no  flow  in  the  plane.  After  it  has  penetrated  the  plane  a 
certain  distance,  it  will  have  a  separated  vortex  flow  pattern. 
The  pattern  shown  in  Figure  2  is  for  a  supercritical  case,  but 
for  a  subcritical  case  the  calculated  crossflow  discrete  vortex 
pattern  would  appear  as  in  Figure  5  for  symmetric  vortices  and  as 
in  Figure  6  for  asymmetric  vortices. 


The  way  in  which  such  vortex  clouds  are  calculated  is  now 
described.  As  the  body  of  revolution  at  some  angle  of  attack 
passes  through  plane  AA,  the  pressure  distribution  around  its 
circumference  in  this  plane  is  calculated,  and  the  lower  stagna¬ 
tion  point  in  the  crossflow  plane  located.  From  this  pressure 
distribution  it  is  possible  to  determine  the  location  of  the 
separation  points  using  the  Stratford  criterion  (assuming 
separation  exists) .  Separation  usually  starts  at  the  top  of  the 
body  for  symmetric  flows.  At  the  separation  points  vortices  are 
introduced  into  the  flow  to  represent  the  vorticity  of  the 
separated  boundary  layer.  This  vorticity  flux  per  unit  time  is 
df/dt  given  by 


dr 

dt 


2  2 


(A) 


Here  is  the  crossflow  velocity  at  the  edge  of  the  boundary 
layer  as  given  by  the  flow  solution  including  any  upstream  vortex 
effects.  The  factor  A  is  an  experimental  factor  equal  to  about 
0.6  to  account  for  vorticity  of  opposite  sign  feeding  into  the 
separation  point  from  above.  A  figure  illustrating  these  points 
is  shown  as  Figure  7.  By  tracking  the  vortices  downstream  and 
applying  the  method  to  a  series  of  crossflow  planes,  the  entire 
vortex  cloud  can  be  calculated.  It  is  noted  that  any  upstream 
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effects  of  the  vortex  cloud  in  influencing  the  pressures  is  not 
included.  This  is  an  approximation  of  the  two-dimensional 
analogy. 

2.4.2  Vortex-Cloud  Method  for  Compressible  Flows.-  The 
incompressible  method  must  be  modified  for  compressibility  in 
several  ways. 

1.  Determine  pressure  distribution  using  Euler  equations. 

2.  Account  for  compressibility  on  separation  position. 

3.  Introduce  compressible  analog  of  incompressible  vortices 
into  the  method. 

With  regard  to  the  first  point,  it  is  of  interest  to  examine 
Figure  8  which  shows  the  flow  in  the  crossflow  planes  of  a  body 
of  revolution  in  supercritical  crossflow  without  vortices.  It  is 
noted  that  crossflow  shocks  going  to  the  body  surface  exist  in  the 
crossflow  plane.  It  is  clear  that  significant  compressibility 
effects  occur  on  the  separation  locations. 

With  regard  to  the  position  of  separation,  it  is  clear  from 
the  foregoing  example  that  new  considerations  have  arisen.  As 
the  boundary  layer  moves  around  the  body  from  the  bottom,  it 
crosses  the  M  =  1  line  and  accelerates  into  the  imbedded  super¬ 
sonic  region.  Whether  it  starts  a  recompression  before  reaching 
the  shock  is  not  clear,  but  at  any  rate  the  separation  mechanism 
may  be  the  shock.  In  this  case  a  different  separation  criterion 
than  that  of  Stratford  may  be  required.  Also  the  flux  of 
vorticity  shed  by  the  separating  boundary  layer  will  probably  be 
different.  Some  theoretical  means  of  handling  this  local  shock- 
wave,  boundary- layer  interaction  will  probably  be  required.  It 
may  possibly  be  obtained  from  existing  shock-wave  boundary-layer 
interaction  methods. 

The  third  question  of  introducing  a  vortex  element  into  the 
flow  to  represent  the  vorticity  flux  of  the  boundary  layer  is  not 
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a  simple  one  aside  from  the  problem  of  superpositioning.  We 
cannot  just  introduce  point  vortices  of  the  incompressible  type 
with  images  inside  the  cylinder.  Since  we  are  calculating  a 
vortical  field  with  a  separation  streamline  emanating  from  the 
surface,  a  crossflow  separation  point  will  be  required.  This 
change  will  introduce  the  desired  flux  into  the  flow  field  at 
the  separation  point. 

In  this  investigation  we  will  consider  the  case  for  which 
the  Mach  number  parallel  to  the  body  axis  is  supersonic  so  that 
we  can  use  a  marching  procedure  for  the  Euler  codes.  Otherwise 
an  extensive  interaction  procedure  will  be  required  using  much 
computer  time.  This  case  will  aiso  avoid  the  shortcomings  of 
the  incompressible  vortex-cloud  method  in  not  accounting  for  the 
upstream  influences  of  the  vortex  cloud. 

It  is  also  clear  that  once  we  can  determine  the  vortical  flow 
using  Euler  equation  methods,  we  can  add  a  wing  and  determine  the 
interference  of  the  body  vortices  on  the  wing.  This  is  the  link 
between  the  previous  work  on  wing-body  interference  without 
forebody  vortices  and  the  present  work. 

A  brief  outline  of  the  general  approach  of  studying  the  body 
vortex  problem  for  steady  supersonic  flows  is  now  given. 

The  steady  Euler  equations  are  written  in  conservation  form 
in  generalized  curvilinear  coordinates.  The  general  coordinate 
system  is  fitted  between  the  body  and  the  outer  bow  shock.  The 
coordinate  system  is  fitted  to  the  body  by  a  conformal  transforma¬ 
tion  which  transforms  the  body  to  a  unit  circle.  To  fit  the 
outer  shock  to  the  coordinate  system,  the  radial  distance  is 
normalized  with  respect  to  the  distance  between  the  outer  shock 
and  the  body.  A  typical  mesh  is  shown  in  Figure  9  for  a  tangent 
ogive  cylinder  body.  The  equations  are  approximated  by 
MacCormack's  second-order  accurate  predictor-corrector  algorithm 
(refs.  13  and  14) .  The  flow  tangency  conditions  at  the  body 
surface  are  satisfied  by  Abbett's  scheme  (ref.  15)  and  the  outer 
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bow-shock  position  by  the  Rankine-Hugoniot  jump  relations.  Any 
internal  shock  waves  and  tangential  discontinuities  or  vortex 
sheets  are  captured  by  the  scheme.  Separated  flow  on  the  body 
will  be  enforced  by  a  Kutta  condition  imposed  at  empirically 
determined  separation  lines. 

In  the  following  section  the  governing  transformed  equations 
will  be  discussed.  Although  these  equations  were  covered  by  the 
previous  final  report  (ref.  3) ,  certain  changes  have  been  made 
and  it  is  worthwhile  to  repeat  these  equations  here  for  com¬ 
pleteness.  The  description  of  the  proper  Kutta  condition  applied 
on  smooth  bodies  follows.  The  final  two  sections  will  cover  the 
numerical  results  and  the  conclusion  and  recommendations.  The 
description  of  the  numerical  scheme,  boundary  and  initial  con¬ 
ditions,  and  body  geometry  will  not  be  covered  in  this  report. 
These  have  been  adequately  covered  by  the  previous  report 
(ref.  3) . 


3.  STEADY  EULER  EQUATIONS 


The  conservation- law  form  of  the  fluid  dynamic  equation  for 
steady,  inviscid,  three-dimensional  compressible  flow  of  an  ideal 
gas  (steady  Euler  equations)  in  Cartesian  form  are 


AAA 


where  E,  F,  and  G  are  defined  as 
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Equation  (1)  represents  the  conservation  of  mass  and 
momentum.  The  pressure  and  density  are  normalized  with  respect 
to  the  stagnation  conditions  and  the  Cartesian  velocity  compo¬ 
nents  (u,v,w)  with  respect  to  the  maximum  adiabatic  velocity 
where  k  =  2Y/(Y-1)  and  Y  is  the  ratio  of  the  specific  heats. 

The  system  of  equations  is  closed  by  the  integrated  form  of  the 
steady  energy  equation  which  in  nondimensional  form  is 

p  =  p(l  -  u2  -  v2  -  w2)  (2) 

For  a  given  free-stream  Mach  number  and  angle  of  attack  a, 
the  remaining  free-stream  variables  in  nondimensional  form  are 
given  by 

,  -Y/(Y-1) 

P«,  *  U  +  [<Y-1)/2JM^} 

,  -1/ (Y-l) 

p„  =  (1  +  ia-D/2]M‘} 


v  =  q  sin  a 
00  **00 

w  =  q  cos  a 

Qo  ^OO 


where 


q  =  (U2  +  V2  +  W2)1/2 
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To  obtain  a  surface-oriented  coordinate  system,  the  system 
(1)  is  transformed  from  the  Cartesian  space  (x,y,z)  into  another 
arbitrary  curvilinear  system  (z,r,<f>)  where  r  =  1  defines  the  body 
surface.  The  general  transformation  is  (the  particular  trans¬ 
formations  considered  are  given  in  section  3  of  reference  3) . 
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(3) 


r  =  r (x,y, z) 

<j>  =  4>  <x,y,  z) 

The  transformed  equations,  obtained  by  Viviand's  procedure 
(described  in  ref.  3) ,  are  in  slightly  rearranged  form 
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where  U,  V,  and  W  are  the  contravariant  velocities  written  with¬ 
out  the  metric  normalization.  The  metric  terms  are  obtained  from 
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where  J  is  the  Jacobian  of  the  transformation  from  the  arbitrary 
curvilinear  space  (z,r ,<f>)  to  the  Cartesian  space  (x,y,z). 

Figure  10.  It  should  be  pointed  out  that  both  coordinate  systems 
are  left-handed,  the  reason  being  that  the  parent  code  (refs.  13 
and  14)  used  for  this  study  was  written  for  a  left-handed  system. 

To  fit  the  outer  bow  shock  wave,  the  outer  mesh  boundary 
must  coincide  with  the  bow  shock.  Since  this  bow  shock  is  a 
variable  three-dimensional  surface  it  is  necessary  to  introduce 
another  transformation  which  normalizes  the  distance  between  the 
body  boundary  and  the  bow  shock  surface.  The  location  of  the  bow 
shock  surface  is  determined  by  the  Rankine-Hugoniot  conditions 
during  the  course  of  the  numerical  computation  described  in  a 
later  section.  At  the  same  time  it  is  desirable  to  have 
arbitrary  clustering  functions  in  the  transformation  so  that 
mesh  points  can  be  concentrated  near  the  body  surface,  wing  tip, 
or  wing-body  juncture  for  increased  resolution  in  areas  of  rapid 
changes  of  the  flow  variables  or  the  previously  mentioned  trans¬ 
formation  metrics. 

The  equations  of  the  independent  variable  transformation 

are 


C  =  z 

x  =  h(Q  (7) 

n  =  f  (<J>) 

where  h  and  f  are  clustering  functions.  The  radial  variable,  £,  is 

5  =  (r  -  rb)/(rs  "  rb> 

where  r^  =  rb(z,<|>),  the  body  surface  radius 

and 

r  =  r  (z,4>),  the  outer  shock  radius 
s  s 


are  to  be  determined  as  part  of  the  numerical  solution  procedure. 
The  system  (1)  may  now  again  be  written  in  strongly  conservative 
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The  functions  h  and  f  are  the  clustering  transformations  in 
the  r  and  <p  directions,  respectively.  The  normalization  between 
the  body  surface  and  the  shock  is  given  by  the  variable  £ . 

The  equation  is  still  in  strongly  conservative  form  unlike 
in  reference  3  where  it  was  thought  that  the  weakly  conservative 
form  is  necessary  to  simplify  the  decoding  of  the  dependent  E 
vector.  The  finite  difference  form  of  equation  (4)  is  integrated 
with  respect  to  the  hyperbolic  coordinate  r,  to  yield  values  of  E. 
The  physical  flow  variables  p,  p,  u,  v,  w  must  be  decoded  from 
the  components  e^  of  E.  Explicit  expressions  are  possible  for 
the  physical  flow  variables  even  if  the  system  is  strongly  con¬ 
servative.  The  decoding  Drocedure  has  been  described  in 


reference  3  and  is  identical  as  required  for  this  equation  set  if 
J  is  replaced  by  JV . 


4.  KUTTA  CONDITION  AND  BOUNDARY  CONDITIONS 

Several  types  of  boundary  conditions  must  be  satisfied  for 
the  problem  considered  in  this  report.  These  include  solid 
surfaces  such  as  the  body  surface,  permeable  surfaces  such  as 
the  fitted  outer  bow  shock  wave,  symmetry  planes,  and  initial 
planes  from  which  the  computation  can  be  started.  Another  type 
of  boundary  considered  is  at  points  on  the  body  where  the  flow 
separates.  At  these  locations  the  flow  tangency  condition  is 
not  adequate.  The  separate  procedure  devised  to  obtain  the  flow 
variables  at  these  points  is  described  below.  The  other  con¬ 
ditions  are  essentially  the  same  as  discussed  in  reference  3  and 
will  not  be  repeated  here. 

In  the  previous  study  (ref.  3)  separated  flows  over  wing- 
body  combinations  were  considered.  The  wings  were  assumed  to  be 
thin  and  planar.  For  these  types  of  wings  with  their  sharp 
leading  edges,  the  separation  line  is  known  a  priori,  namely  at 
the  sharp  leading  edges.  Furthermore  it  was  shown  previously 
that  at  the  sharp  leading  edges  the  transformation  metrics  and 
Jacobian  vanish.  Thus  arbitrary  values  of  the  flow  variables 
will  satisfy  the  Euler  equation  since  these  equations  are 
indeterminate  at  the  leading  edges.  However,  for  the  tangent 
ogive  cylinders  considered  in  this  report,  the  separation  line 
is  not  known  a  priori.  Neither  do  the  Euler  equations  become 
indeterminate  anywhere  on  the  body.  Thus  a  different  type  of 
Kutta  condition  must  be  prescribed  for  these  types  of  bodies. 

The  specification  of  the  Kutta  condition  to  enforce  separated 
flow  is  straightforward.  The  five  flow  variables  must  be 
specified.  They  are  the  three  Cartesian  velocity  components. 
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u,  v,  w  and  the  pressure  and  density.  The  values  of  these  five 
variables  are  not  arbitrary  as  in  the  previous  study  (ref.  3) 
since  the  transformed  equations  are  not  indeterminate.  This 
means  that  the  flow  variables  must  satisfy  the  continuity, 
momentum,  and  energy  equations.  Since  we  want  to  force  separa¬ 
tion,  the  flow  variables  must  satisfy  the  conditions  for  a 
tangential  discontinuity.  The  conditions  for  a  tangential  dis¬ 
continuity  surface  or  vortex  sheet  are 

a)  p  is  continuous  across  the  sheet 

b)  the  velocity  normal  to  the  sheet  is  zero 

c)  total  enthalpy  is  constant 

d)  the  jump  in  the  velocity  tangential  to  the  sheet  is 
arbitrary 

e)  the  jump  in  the  density  p  across  the  sheet  is 
arbitrary 

The  jump  relations  state  that  the  jumps  in  density  and  tangential 
velocity  are  arbitrary.  This  means  that  the  jumps  are  not 
determined  by  the  local  conditions,  but  by  global  conditions. 

The  global  conditions  are  determined  by  the  governing  equations 
and  boundary  conditions  and  not  the  jump  relations. 

The  vortex  sheet  coming  off  the  body  is  shown  in  Figure  11. 
Two  geometric  parameters  determine  the  location  and  orientation 
of  this  sheet.  They  are  4>  and  <J>  for  the  location  and  orienta- 
tion,  respectively,  as  shown  in  Figure  11.  On  each  side  of  this 
vortex  sheet,  separate  values  of  p,  p,  u,  v,  w  must  be  prescribed. 
However,  since  this  sheet  is  to  be  captured  by  the  finite  differ¬ 
ence  scheme  it  is  more  appropriate  to  prescribe  some  mean  values 
of  p,  o,  u,  v,  w  within  the  sheet.  These  mean  values  will  be 
some  weighted  average  of  the  values  on  either  side  of  the  sheet. 
Thus  we  can  prescribe  some  mean  velocity  vector  within  the  sheet 


20 


by  qQ  (Fig.  11) .  This  is  the  mean  tangential  velocity  vector 
oriented  by  angle  <f>  from  the  axial  direction  (z-coordinate)  . 

The  conditions  for  separation  at  a  point  on  the  body  4>  =  <p 
are  implemented  as  follows,  with  reference  to  Figure  12: 

a)  pQ  =  0.5(p1  +  p_1) 

b)  pQ  =  0.5(p1  +  p^) 

c)  total  enthalpy  is  constant 

where  the  subscript  "o"  designates  the  separation  point,  "1"  the 
first  point  on  the  body  on  the  leeside  of  the  vortex  sheet,  and 
"-1"  the  first  point  on  the  body  on  the  windward  side.  In  other 
words,  the  pressure  and  density  at  the  separation  point  are 
obtained  by  averaging  over  the  nearest  surface  neighboring  values. 
These  three  conditions  are  sufficient  to  determine  the  mean 
tangential  velocity  qQ  in  the  vortex  sheet  as  shown  in  Figures  11 
and  12.  To  obtain  the  Cartesian  velocity  components  uq,  v  ,  and 
w  ,  the  three  geometric  parameters  <j>  ,  <J>  ,  and  <p  are  required. 
Since  the  Euler  equations  do  not  account  for  the  viscous  effects, 
which  control  the  separation,  additional  information  is  required 
to  determine  the  three  geometric  parameters.  These  parameters 
will  be  determined  from  empirical  data  for  the  present  study. 

The  three  velocity  components  at  the  separation  point  are  given 
by 
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The  above  procedure  has  been  used  successfully  in  computing 
separated  flow  over  smooth  bodies.  However,  it  may  be  possible 
to  improve  the  procedure.  For  example,  condition  (a)  simply 
averages  the  pressure  over  its  nearest  neighbors.  Some  experi¬ 
mental  data  seem  to  indicate,  though,  that  the  pressure  reaches 
a  local  extremum  at  the  crossflow  separation  points.  This  can  be 


i 


21 


seen  in  Figure  13  (Peake  and  Tobak,  ref.  16),  which  shows  the 
location  of  the  separation  points  by  an  oil  film  and  the 
corresponding  pressure  in  the  crossflow  plane.  A  simple  pressure 
averaging  would  thus  underpredict  the  pressure  at  the  separation 
point.  A  more  accurate  method  would  be  to  extrapolate  the  pres¬ 
sures  to  the  separation  point  from  each  side  and  then  average. 
This  has  not  been  tried  in  this  study. 

It  would  also  be  useful  to  obtain  the  three  geometric 
parameters  from  some  sort  of  separation  criteria  such  as  the 
Stratford  criteria  which  works  well  for  subsonic  flows.  If  a 
compressible  separation  criteria  can  be  established,  then  the 
Euler  code  can  be  used  to  determine  separated  flows  independently 
of  any  particular  set  of  empirical  data. 


5.  NUMERICAL  RESULTS 

In  this  section  some  of  the  numerical  results  obtained  by 
the  Euler  code  with  a  Kutta  condition  will  be  covered.  Two  cases 
are  discussed.  Firstly,  a  two-caliber  tangent-ogive  cylinder  at 
a  free-stream  Mach  number  of  3.0  and  angle  of  attack  (a)  of  15° 
is  presented.  This  will  give  a  direct  comparison  to  the  experi¬ 
mental  data  of  Oberkampf  (ref.  6)  and  to  the  vortex-cloud  method 
of  Mendenhall  (ref.  17) .  Secondly,  a  five-caliber  tangent-ogive 
cylinder  at  M^  =  2.9  and  a  =  20°  is  computed.  This  case  will 
allow  comparison  with  results  obtained  by  a  parabolized  Navier- 
Stokes  code  by  Pulliam  (ref.  18) .  A  parametric  study  of  the 
geometric  parameters  required  for  the  Kutta  condition  is  also 
carried  out.  These  results  will  be  presented  first. 

5.1  Geometric  Parametric  Study 

A  parametric  study  of  the  three  geometric  parameters  <J>  ,  4>  , 

s  c 

and  6  is  required  to  complete  the  specification  of  the  Kutta 

cl 
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condition.  This  will  determine  the  relative  importance  of  each 
of  the  parameters  in  controlling  the  size  of  the  separated  flow 
region  and  the  amount  of  vorticity  shed  into  this  region.  This 
study  was  carried  out  on  a  three-caliber  tangent-ogive  nose  with 
a  cylindrical  afterbody  at  =  3.0  and  a  =  20°.  A  simple  way 
of  determining  the  relative  effects  of  the  parameters  is  to 
perturb  each  of  the  parameters  around  a  base  combination  of 
<b  =  90°,  <b  =  30°,  and  <4  =  10°.  The  base  combination  values 

SC  a 

were  estimated  from  experimental  data  such  as  those  shown  in 
Figure  13  (obtained  from  ref.  16) . 

The  results  indicate  that  the  most  influential  parameter  is 

the  separation  angle,  <t>s.  This  parameter  determines  the  location 

of  the  separation  vortex  sheet  on  the  body.  As  shown  in  Figure 

14(a),  the  separation  angle  has  a  significant  effect  on  both  the 

total  amount  of  vorticity  shed  into  the  separated  region  and  on 

the  height  of  this  region  as  measured  on  the  leeward  plane  of 

symmetry.  The  crossflow  angle,  c ,  seems  to  have  very  little 

effect  on  the  rate  of  vorticity  shedding,  Figure  14(b).  Both  the 

vertical  height  of  the  leeward  plane  stagnation  point  and  the 

total  circulation  are  essentially  independent  of  this  parameter. 

The  axial  angle  <t  has  a  rather  moderate  effect  on  the  vorticity 

« 

shedding  rate  and  a  somewhat  stronger  effect  on  the  height; 

Figure  14  (c) . 

The  results  indicate  that  the  major  influential  geometric 

parameter  is  the  separation  angle  <J)g.  The  axial  angle  parameter 

is  of  only  moderate  importance  and  the  crossflow  angle  <b  seems 

c 

to  have  virtually  no  influence.  This  means  that  if  the  param¬ 
eters  are  to  be  determined  empirically,  the  most  accuracy  required 
is  on  4>s  and  rough  estimates  are  sufficient  for  the  remaining  two. 
This  is  fortuitous  since  the  separation  angle  can  be  accurately 
determined  from  oil  flow  or  vapor  screen  data.  However,  it  is 
more  difficult  to  determine  the  axial  angle  <J> 

a 

angle  <f>  . 

c 


and  the  crossflow 


5.2  Pressure  Interpolation 


As  mentioned  in  section  4  the  pressure  at  the  separation 
line  is  determined  by  linear  interpolation  from  the  nearest 
neighboring  points  on  the  body.  It  is  not  clear  that  this  is 
the  best  strategy.  Although  the  pressure  is  constant  across  a 
vortex  sheet,  there  may  be  strong  pressure  gradient  leading  up 
to  the  separation  line  as  indicated  by  Peake  and  Tobak  (ref.  16) 
and  shown  in  Figure  13.  If  the  mesh  spacing  is  fine  enough  near 
the  separation  line,  then  linear  interpolation  should  be  adequate. 
Some  pressure  data  are  available  from  reference  7  with  which 
direct  comparisons  can  be  made  for  the  numerical  calculation 
using  the  base  values  of  the  geometric  parameters  of  the  previous 
section.  These  comparisons  are  given  in  Figures  15a  through  15c 
for  three  different  axial  stations.  Shown  are  the  circumf erential 
pressure  distributions  for  a  3.0-caliber  tangent-ogive  nose  with 
a  cylindrical  afterbody  at  a  =  20°  and  =  2.96  for  the  experi¬ 
mental  data  and  M  =3.00  for  the  numerical  results.  The 
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agreement  is  reasonable  except  at  the  windward  crossflow  stagna¬ 
tion  point  (<t>  =  0°)  and  over  the  separated  region.  There  are 
probably  two  reasons  for  this  discrepancy.  The  numerical  results 
have  the  Kutta  condition  imposed  at  <p  =  90°  which  is  the  approxi¬ 
mately  asymptotic  value  of  the  experimental  separation  line. 

The  experimental  separation  line  varies  as  shown  in  Figure  16 
from  <(>  =  160°  at  z/D  =  2.0  to  4>g  =  90°  at  z/D  =  80.  This 
separation  line  is  determined  from  the  oil  flow  data  shown  in 
Figure  17.  The  asymptotic  value  at  z/D  2  8  is  determined  from 
similar  data  from  Peake  and  Tobak  (ref.  16) .  Thus  the  separation 
occurs  at  different  circumferential  locations  and  the  discrep¬ 
ancies  are  to  be  expected.  The  other  reason  is  probably  the  lack 
of  accuracy  in  determining  the  axial  angle,  $  .  While  this 
parameter  has  only  a  moderate  effect  on  the  overall  flow  separa¬ 
tion,  i.e.  size  and  vortex  shedding  rate,  it  seems  to  have  a  much 
stronger  effect  on  the  local  pressure  distribution.  This  can  be 
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seen  in  Figure  15  as  a  large  oscillation  in  the  pressure  near  the 
imposed  separation  point  (4>  =  90°)  for  an  axial  angle  4>a  =  10°. 

This  case  was  rerun  with  a  variable  separation  line  given  in 
Figure  16  from  the  oil  flow  data  of  Figure  17.  With  the  axial 
angle  varying  from  0°  at  z/D  =  2.0  to  only  5°  at  z/D  =  8.0  the 
results  shown  in  Figures  18a  and  b  were  obtained.  The  large 
oscillations  are  nearly  absent  but  agreement  with  the  experi¬ 
mental  data  is  not  notably  improved  in  the  separated  region. 

There  may  be  two  or  more  separation  lines.  Unfortunately  the 
experimental  pressure  data  resolution  is  rather  coarse  (every 
22.5°  circumferentially)  for  this  case.  Thus  it  is  not  possible 
to  discern  the  additional  separation  lines.  The  oil  flow  data 
are  also  not  clear  enough  to  determine  these.  Although  the  pres¬ 
sure  comparisons  may  be  suspect,  the  overall  flow  field  prediction 
is  good  such  as  the  location  of  the  outer  bow  shock  and  the  over¬ 
all  size  and  shape  of  the  separated  region.  This  may  be  seen  in 
Figure  19,  which  shows  the  numerical  results  superimposed  on  the 
experimental  data.  It  is  somewhat  difficult  to  discern  the  size 
of  separated  zone  from  the  schlieren  and  vapor-screen  photographs. 
For  this  reason  another  set  of  runs  was  made  for  a  2.0-caliber 
tangent-ogive  cylinder  for  which  detailed  flow  field  data  are 
available . 

5.3  Numerical  Results 

Other  numerical  results  have  been  obtained  for  two  tangent 
ogive  circular  bodies  at  two  angles  of  attack  and  free-stream 
Mach  numbers.  The  first  case  for  which  there  exist  extensive 
experimental  flow  field  data  within  the  separated  region  is  for 
a  2.0  caliber  tangent-ogive  nose  and  cylinder  combination  at 
Mach  3.0  and  angle  of  attack  of  15°.  The  second  case  is  a  5.0 
caliber  tangent-ogive  nose  and  cylinder  combination  at  M^  =  2.9 
and  a  =  20°.  For  this  latter  case  an  independent  numerical 
calculation  based  on  the  parabolized  Navier-Stokes  equation 
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(ref.  18)  is  available.  For  both  of  these  cases  the  Kutta  con¬ 
dition  is  started  at  z/D  =  2.0  and  the  separation  angle  $  is 
kept  constant  at  90°.  The  vortex  orientation  angle,  <J>  ,  and  the 
axial  angle,  <p  ,  are  30°  and  10°,  respectively. 


5.3.1  Two  Caliber  Tangent  Ogive  Cylinder.-  For  this  con¬ 
figuration  direct  comparison  may  be  made  with  the  extensive 
experimental  data  of  Oberkampf  (refs.  4,  5,  and  6).  These  data 
cover  only  the  flow  field;  surface  data  such  as  circumferential 
pressure  distribution  and  oil  flow  data  are  not  available.  The 
first  numerical  result  presented  are  the  crossflow  velocity 
vectors  at  a  station  thirteen  calibers  downstream  of  the  nose 
tip.  Figure  20.  Only  the  leeward  quadrant  is  shown.  As  can  be 
seen  from  the  figure,  separation  occurs  at  $  =  90°  and  the  entire 
lee  side  of  the  body  is  in  the  reversed  flow  field.  No  allowance 
was  made  for  a  possible  secondary  separation.  Other  experimental 
data  (Peake  and  Tobak,  ref.  16)  indicate  that  for  this  Mach 
number,  angle  of  attack  and  configuration,  secondary  separation 


occurs.  Also  indicated  in  this  figure  is  the  approximate  loca¬ 
tion  of  the  vortex  center  and  the  height  (by  ystag)  of  the 
separated  flow  region.  The  heavy  dotted  line  outlines  the  area 
covered  by  Oberkampf' s  data  in  Figure  21.  The  experimental 
vortex  center  is  somewhat  smeared  out  but  overlaps  the  predicted 


location.  The  extent  of  the  separated  region  is  predicted  within 


3  percent  of  the  body  diameters.  Unfortunately  the  surface  data 
are  not  available  for  determining  the  primary  or  secondary  separa¬ 
tion  point.  The  axial  growth  of  the  circulation  around  the  contour 
shown  in  Figure  20  by  the  heavy  dotted  line,  plane  of  symmetry, 
body,  and  <p  =  90°  is  shown  in  Figure  22.  The  agreement  with 
Oberkampf' s  data  is  good.  Also  shown  are  the  results  obtained 
by  the  vortex-cloud  theory  of  Mendenhall  (ref.  17).  It  con¬ 
sistently  underpredicts  the  circulation  for  this  Mach  number. 


The  crossflow  particle  paths  of  the  numerical  prediction  are 
depicted  in  Figure  23.  This  plot  shows  the  flow  field  as  would 


be  seen  by  a  vapor-screen  visualization.  There  is  a  spurious 
vortex  just  behind  the  separation  line  which  is  due  to  the  large 
surface  pressure  oscillation.  The  main  vortex  is  typically  of 
those  observed  experimentally,  i.e.  elongated  and  elliptical  in 
shape.  It  is  possible  that  the  spurious  vortex  indicate  that 
the  numerical  results  should  have  allowed  for  a  secondary 
separation. 

While  Oberkampf  did  not  provide  any  vapor-screen  visualiza¬ 
tions,  he  did  provide  detailed  data  from  which  velocity  profile 
through  the  separated  region  can  be  plotted.  The  first  series, 
Figure  24,  are  the  vertical  velocity  profiles  along  the  lee 
plane  of  symmetry  for  several  z-stations.  The  size  of  the 
separated  region  (i.e.  Ystag  =  y  where  v/V^  =  0)  is  predicted 
fairly  well,  but  the  peak  of  the  maximum  vertical  velocity  com¬ 
ponent  is  underpredicted.  It  is  not  a  radial  mesh  resolution 
problem  since  an  increase  of  the  number  of  radial  points  (between 
the  body  and  bow  shock  wave)  from  30  to  40  did  not  improve  the 
results.  However,  it  could  be  that  not  enough  circumferential 
points  are  in  the  separated  region.  This  seems  to  be  borne  out 
by  the  following  three  sets  of  figures.  These  figures  show  the 
vertical  velocity  profiles  along  horizontal  lines.  Figure  25  is 
for  the  horizontal  line  at  y/D  =  0.3.  Figures  26  and  27  are  for 
y/D  =  0.6  and  0.9,  respectively.  In  all  these  curves  the  pro¬ 
files  seem  to  be  excessively  smeared.  This  is  indicative  of  not 
enough  mesh  resolution,  which  has  the  tendency  to  reduce  the 
sharp  velocity  gradients  in  the  recirculating  region. 

5.3.2  Five  Caliber  Tangent-Ogive  Cylinder.-  The  second  case 
considered  in  this  section  is  a  5  caliber  tangent-ogive  nose  and 
cylinder  combination  at  =  2.9  and  a  =  20°.  For  this  case  the 
numerical  computations  of  Pulliam  (ref.  18)  based  on  the 
parabolized  Navier-Stokes  (PNS)  equation  are  available.  The  PNS 
code  is  based  on  the  steady  Navier-Stokes  equations.  For  a 
supersonic  axial  velocity  the  equations  are  hyperbolic  in  that 
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direction  and  a  marching  procedure  is  employed.  In  the  boundary- 
layer  region  where  the  velocity  becomes  subsonic  special  tech¬ 
niques  are  employed  to  maintain  the  hyperbolic  form  of  the  Navier- 
Stokes  equations  over  the  entire  flow  field.  The  thin  layer 
model  is  also  employed,  which  means  that  only  the  viscous  dif¬ 
fusion  in  the  direction  normal  to  the  body  surface  is  retained. 
Diffusion  terms  in  the  streamwise  and  crossflow  direction  are 
neglected. 

The  first  set  of  results  (Figs.  28(a)  and  28(b))  shows  the 
crossflow  velocity  vectors  from  the  Euler/Kutta  code  and  the  PNS 
code  computations.  As  expected  there  are  large  differences  near 
the  surface  due  to  the  viscous  effects.  What  is  not  expected 
however  are  the  large  differences  between  the  size  and  center 
location  of  the  separated  region.  The  Euler/Kutta  code  predicts 
a  vortex  center  at  about  y/D  *  1.2  and  the  PNS  code  at  about 
y/D  ~  0.7.  Similarly  the  heights  of  the  recirculating  region 
are  Ystag  a  2.0D  and  1 . ID  for  the  Euler/Kutta  and  PNS  codes, 
respectively.  These  results  seem  to  cast  some  serious  doubts 
on  the  Euler/Kutta  procedure.  However,  comparing  the  two 
numerical  results  with  extensive  experimental  data  (ref.  12) 
reveals  that  the  Euler/Kutta  results  falls  within  the  expected 
range  of  values  determined  experimentally.  Figure  29,  for  the 
height  of  the  recirculation  region  at  z/D  =  10.  The  PNS  code 
result,  on  the  other  hand,  completely  underpredicts  the  height 
of  the  recirculating  region. 

The  main  cause  for  this  discrepancy  is  probably  the  lack  of 
resolution  for  the  PNS  code  results  in  the  region  of  the  separated 
flow  outside  the  boundary  layer.  Although  the  resolution  is  not 
that  much  better  for  the  Euler/Kutta  results,  lack  of  resolution 
is  much  more  important  for  the  PNS  code.  This  code  is  based  on 
an  implicit  scheme  which  has  a  much  higher  dissipation  inherent  in 
the  algorithm  than  the  explicit  MacCormack's  scheme  used  for  the 
Euler/Kutta  code.  Thus  the  implicit  scheme  dissipation  is  much 
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more  mesh  size  dependent  than  the  MacCormack  scheme.  This  may 
be  possibly  another  advantage  in  the  Euler/Kutta  procedure. 

Namely,  reasonably  good  results  are  obtained  with  a  fairly 
coarse  mesh. 

A  comparison  of  the  execution  time  for  the  two  codes  can  be 
made.  The  Euler/Kutta  code  requires  about  1000  seconds  on  the 
CDC  7600  computer  to  advance  the  solution  to  a  station  16 
calibers  downstream  of  the  nose  tip  for  a  mesh  of  30  radial 
and  36  circumferential  points.  The  PNS  code  requires  about  1100 
seconds  to  advance  the  solution  to  the  same  streamwise  stations 
for  a  mesh  of  40  radial  and  24  circumferential  points.  It  is 
clear,  however,  that  the  circumferential  resolution  for  the  PNS 
code  is  inadequate.  If  the  comparison  is  made  on  the  basis  of 
the  same  circumferential  resolution,  then  the  PNS  code  requires 
about  1600-1700  seconds.  This  assumes  that  the  radial  resolution 
of  40  points  is  adequate.  Thus  the  PNS  code  requires  about  60 
to  70  percent  more  computer  time  than  the  Euler/Kutta  procedure, 
even  though  the  implicit  algorithm  is  more  efficient.  That  is, 
it  has  no  stability  limitations  on  the  marching  step  size. 

6.  CONCLUSIONS  AND  RECOMMENDATIONS 

This  report  addressed  the  problem  of  computing  the  steady 
inviscid  rotational  supersonic  flows  about  tangent  ogive  cylinders 
At  high  angles  of  attack  the  flow  separates  into  a  pair  of  counter 
rotating  vortices  on  the  lee  side  of  the  cylinder.  These  vortices 
have  an  adverse  effect  on  the  lift  of  any  wings  or  fins  attached 
to  the  cylinder.  Previous  studies  of  the  body  vortices  include 
the  vortex-cloud  method  of  Mendenhall  (ref.  12)  and  numerical 
calculations  using  the  parabolized  Navier-Stokes  equation  by 
Pulliam  (ref.  18) .  The  vortex-cloud  method  is  based  on 
Stratford's  separation  criterion  which  is  not  applicable  if  the 
crossflow  Mach  number  is  greater  than  0.6.  The  numerical 
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calculations  of  Pulliam  are  rather  lengthy  even  though  they  are 
based  on  the  parabolized  version  of  the  more  complete  Navier- 
Stokes  equations. 

We  have  calculated  separated  flows  with  symmetric  vortices 
on  tangent  ogive  cylinders.  The  inviscid  steady  Euler  equations 
have  been  used  as  the  basic  governing  equations  as  opposed  to  the 
more  complicated  Navier-Stokes  equations  to  obtain  a  more  efficient 
procedure.  The  viscous  effects,  important  only  near  the  separa¬ 
tion  lines,  are  simulated  by  a  Kutta  condition.  The  rest  of  the 
flow  field  is  essentially  controlled  by  the  inviscid  equations. 

The  steady  Euler  equations  were  written  in  conservation  form 
in  generalized  curvilinear  coordinates.  The  general  coordinate 
system  was  fitted  to  a  cylindrical  body  by  a  conformal  transforma¬ 
tion  between  the  body  and  a  unit  circle.  The  outer  radial  mesh 
points  were  adjusted  to  the  bow  shock  wave  by  normalizing  the 
radial  coordinate  between  the  body  and  the  shock.  Although  only 
circular  bodies  were  considered  in  this  report,  the  conformal 
transformations  will  facilitate  the  application  of  this  procedure 
to  noncircular  bodies  and  wing-body  combinations  later.  The 
equations  were  approximated  by  MacCormack ' s  second-order  accurate 
predictor-corrector  algorithm  (ref.  14)  with  special  stabilizing 
terms  appended.  The  flow  tangency  conditions  at  the  body 
surface  were  satisfied  by  Abbett’s  scheme  (ref.  15)  and  the 
outer  bow  shock  position  by  the  Rankine-Hugoniot  jump  relations. 
Any  internal  shock  waves  or  tangential  discontinuities  were 
captured  by  the  scheme . 

The  Euler  equations  and  Kutta  condition  approach  was  used 
in  a  previous  study  of  planar  delta  wings  and  wing-body  combina¬ 
tions  (ref.  3) .  For  these  configurations  with  sharp  leading 
edges  the  separation  line  was  known  a  priori.  For  the  tangent 
ogive  cylinder  considered  in  this  report,  the  separation  line 
could  not  be  determined  from  inviscid  theory.  For  this  reason. 
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empirical  data  were  required  to  determine  the  separation  lines. 

A  parametric  study  of  the  various  parameters  that  affect  the  flow 
separation  lines  off  the  body  indicated  that  the  overall  flow 
field  is  relatively  insensitive  to  the  fine  details  near  the 
separation  point.  The  most  important  parameter  was  the  location 
of  the  separation  line  on  the  body.  The  orientation  of  the 
tangential  discontinuity  sheet  relative  to  the  body  is  of  only 
secondary  importance  in  controlling  the  overall  flow  field.  The 
question  of  determining  the  pressure  and  density  at  the  body 
separation  points  is  still  open.  But,  nevertheless,  good  results 
were  obtained  with  simple  linear  interpolation  and  with  a  fairly 
coarse  mesh.  A  finer  mesh  produced  better  overall  results  and  a 
smoother  pressure  distribution  on  the  body.  It  turned  out  that 
specifying  a  variable  separation  line  on  the  nose  region  did  not 
greatly  influence  the  results  further  downstream  as  long  as  the 
proper  asymptotic  value  of  the  separation  angle  was  finally 
reached.  Specifying  the  asymptotic  value  along  the  entire  length 
produced  similar  results.  The  variable  separation  line  produced 
only  minor  coding  problems.  A  second  Kutta  condition  could  be 
applied  to  obtain  the  secondary  separated  flow.  However,  it  is 
not  known  how  important  this  secondary  separation  may  be.  It  is 
suspected  that  it  will  give  better  surface  pressure  results  in 
the  separated  flow  region  and  have  very  little  influence  else¬ 
where  . 

It  was  also  shown  that  the  Euler/Kutta  codes  gave  results 
which  agreed  well  with  experimental  data.  These  results  were 
obtained  for  much  lower  run  times  than  those  from  the  PNS  code. 

It  is  anticipated  that  these  costs  could  be  further  reduced  (by 
a  factor  of  two  or  three)  if  the  explicit  algorithm  used  in  the 
code  is  replaced  by  the  more  efficient  implicit  codes  used  for 
the  PNS  code . 

We  have  shown  that  the  Euler/Kutta  procedure  works  well  for 
smooth  bodies  such  as  cylindrical  bodies  with  tangent-ogive 
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noses  considered  in  this  report  as  long  as  empirical  data  are 
available  to  specify  a  priori  the  separation  line.  The  applica¬ 
tion  of  this  procedure  for  sharp-edge  bodies  such  as  wings  has 
been  demonstrated  in  the  previous  year's  work.  Therefore  we  are 
now  in  a  position  to  compute  the  interaction  problem  of  a 
complete  missile  configuration.  This  problem  can  include  the 
interaction  of  the  forebody  vortices  with  the  wings  and  the  wing 
vortices  as  well  as  wing-tail  interference. 

However,  before  such  a  capability  to  solve  a  complete 
missile  configuration  is  established,  it  is  recommended  that  the 
following  work  be  carried  out. 

1.  The  determination  of  the  pressure  and  density  at  the 
separation  line  on  the  body  needs  further  clarification. 

2.  The  implicit  algorithm  should  be  adapted  for  the  Euler/ 
Kutta  code.  This  will  further  reduce  the  computation 
time  by  a  factor  of  two  or  three. 

3.  The  effects  of  including  secondary  separation  points  on 
the  forebody  vortices  should  be  assessed. 

4.  Finally,  a  separation  criterion  for  supercritical  cross- 
flow  should  be  established  from  empirical  data. 
Incorporation  of  this  separation  criterion  will  make 
the  Euler/Kutta  code  independent  of  any  specific  set 

of  empirical  data  and  thus  much  more  useful. 
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(a)  Symmetric  separation 


(b)  Asymmetric  separation 


1.-  Symmetric  and  asymmetric  vortex  formation 
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Figure  3.-  Estimated  height  of  top  of  recirculation  region 

at  X/D  =  10  for  tangent-ogive  cylinder  (ref.  2). 


Figure  5.-  Calculated  symmetric  vortex  cloud  in  crossflow 
plane  of  L/D  =  7  ogive-cylinder. 


Figure  6.-  Calculated  asymmetric  vortex  cloud  in  crossflow 
plane  of  L/D  =  7  ogive-cylinder. 
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Vortical  singularity 


Figure  8.-  Inviscid  supersonic  crossflow  about 
body  at  high  angle  of  attack. 
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Transformed  Space 


Figure  12.-  Detail  of  separation  point  in 
the  cross  flow  plane. 
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Normalized  circulation,  r/ (ttdv 


Figure  14.-  Effect  of  the  three  geometric  parameters  on  the 
vorticity  and  size  of  the  separated  flow  region  of  a 
3.0  caliber  tangent  ogive  cylinder 
at  z  =  11  calibers. 


height  of  leeward  recirculating  region 


Normalized  circulation,  r/(TDV( 


(b)  Crossflow  angle,  <f> 

C 

Figure  14.-  Continued. 
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height  of  leeward  recirculating  region 


Normalized  circulation,  P/(TOV 


(c)  Axial  angle,  $ 
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Figure  14.-  Concluded. 


50 


height  of  leeward  recirculating  region 


(b)  z/D  =  5.0 
Figure  15.-  Continued. 
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.  20 


Landrum's  data  (Ref.  7) 

L  /D  =  3/9  L/D  =6.67 


(a)  z/D  =  4.5 


Figure  18.-  Circumferential  pressure  distribution  on 
a  tangent-ogive  cylinder  at  a  =  20°. 
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tangent  ogive  nose 


cloud  theory,  and  experimental  data. 
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Figure  23.-  Cross  flow  particles  paths  as  determined 
by  Euler  code  with  Kutta  condition. 
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Figure  24.-  Vertical  velocity  component  on  lee  plane 
of  symmetry  (cj>  =  180°)  . 


Figure  24.-  Concluded 
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Figure  25.-  Continued 


Figure  25.-  Concluded 


Figure  26.-  Vertical  velocity  component  vs.  x/D  at  y/D 


Figure  26.-  Concluded 


Vertical  velocity  component 


Figure  27.-  Continued 
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(C)  z/D  =  13.0 
Figure  27.-  Concluded 
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Figure  29.-  Estimated  height  of  top  of  recirculation  region 
at  z/D  =  10  for  tangent-ogive  cylinder. 
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